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METHODS OF DETERMINING THE ORBITS OF 
SPECTROSCOPIC BINARIES. 



By Heber D. Curtis. 



With the constantly increasing number of discoveries in this 
field of spectrograph^ research, the subject of spectroscopic 
binaries as a whole has become an exceedingly important one 
in its bearing on theories of stellar evolution. Thirty, or even 
twenty, years ago it would have seemed beyond the bounds 
of reasonable conjecture that at least one in every seven or 
eight of the brighter stars should prove to be a double star 
whose components are impossible of separation in any tele- 
scope existing or whose construction is conceivable. It is 
unnecessary to speak of the importance of the results already 
secured in relation to theories of cosmogony. As one con- 
clusion only, we are led to think that our Sun, with its retinue 
of small planets, may possibly be a rather exceptional case in 
stellar development, the more common case being the fission 
of a star into two or more masses of the same order of 
magnitude. 

The number of these systems is increasing to-day at a rapid 
rate. Including those whose binary character is under inves- 
tigation, but not sufficiently substantiated for publication, the 
number is now about two hundred; and, with the carrying of 
spectrographic surveys down to fainter stars, we may need 
four figures to express the total in the not distant future. 

It seems improbable at present that the labor of "keeping 
track" of all these systems will ever present the difficulties 
which make the asteroid group rather an embarrassing source 
of astronomical wealth, though we cannot, of course, premise 
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just what interesting developments may be gained as to the 
inner mechanism of these systems from a comparison of their 
orbits of to-day with their forms one hundred years hence. 
It is not too much to suppose that the student of celestial 
mechanics may then find in this field "laboratory demonstra- 
tions" of the problem of three or more bodies, while possible 
changes in the periods and eccentricities could be expected to 
give valuable results to the worker in theories of tidal evolution 
of planets and satellites. 

In the sphere of double stars we now have an unbroken 
continuum, from the spectroscopic binary of a few hours' 
rotation period to the visual double star with a period so long 
that centuries will be necessary for an accurate determination 
of the orbital elements. Some five or six stars are known to 
be both spectroscopically and telescopically double, and these 
serve to bind together the two fields. Visual double-star 
observers to-day are agreed that the greatest present progress 
in their branch of this wide field is to be derived from the 
systematic measurement of the systems of shorter period, from 
five to fifty years. It is in just this point of shortness of period 
that the spectroscopic field has the advantage in the point of 
evolutional developments. To take an extreme case: In the 
interesting quadruple system of Castor the published spec- 
trographic observations of a u the fainter component, extend 
over eleven years and cover no less than fourteen hundred 
complete revolutions, while for the visual system the telescopic 
measures extending over one hundred years cover so short an 
arc that the visual elements must remain uncertain for a long 
time to come. 

The actual determination of the elements of a spectroscopic 
binary is- comparatively easy. From twenty to thirty plates, 
if well distributed over the period of revolution, are ample, in 
most cases, to secure satisfactory values of the elements of the 
system. The securing of these plates is the time-consuming 
factor in the problem, however, which perhaps explains why, 
for about seventy-five per cent of known spectroscopic binaries, 
we are still without accruate elements. While a skilful 
observer of visual double stars can make and reduce a com- 
plete single set of measures in fifteen to twenty minutes, 
including in this the time taken in passing from one star of 
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his programme to another, from half an hour to two or more 
hours is required for the exposure time alone of a single 
spectrographic plate. Developing, measuring, and reducing, 
unless but few lines are available in the spectrum, will take 
several hours more. From a consideration of these points it 
may easily be seen that the determination of the orbits of 
three or four hundred spectroscopic binaries would require in 
the neighborhood of ten thousand plates, — easily twenty years' 
work for a single observatory, however favorably located as to 
climatic conditions. 

Doubtless for some years to come the individual observa- 
tories possessing powerful spectrographs for radial-velocity 
determinations will endeavor to compute the elements of the 
binaries discovered by themselves. But it seems probable that 
the time must come when there will be systematic co-operation, 
to some extent at least, in order to secure greater rapidity in 
utilizing the data at present available and to be accumulated, 
and to prevent unnecessary duplication of work. 1 

As no recent treatment of the various methods for the 
determination of the elements of a spectroscopic binary exists 
in English, it is the purpose of this paper merely to collect in 
accessible form some of the more approved methods of solving 
the problem, with a few hints as to the processes involved, but 
without pretending to advance any new matter on the subject. 

The more important papers on the subject may be referred 
to as follows: — 

Rambaut, M. N., Vol. LI, pp. 316 ff. 
Wilsing, A. N., 134, 90, 1893. 
Russell, A ph. J., 15, 252, 1902. 
Lehmann-Filhes, A. N., 136, 17, 1894. 

SCHWARZSCHILD, A. N ., 152, 65, I9OO. 

Bauschinger, "Bahnbestimmung der Himmelskorper," 
pp. 647 ff. 

ZURHELLEN, A. N ., 173, 353, I9O7. 
ZlJRHELLEN, A. N ., 175, 245, I907. 



1 For several years the Lick Observatory has been keeping records of work in 
this field, — a task necessitated by the need for systematizing its own efforts, — and 
as a result it was found wise to collect the material available in a "First Catalogue 
of Spectroscopic Binaries" (Campbell and Curtis, Lick Observatory Bulletin, 
No. 79, 1905). From time to time, as the need arises, successive catalogues will 
be issued. — W. W. C. 



136 Publications of the 

Of these methods, that of Rambaut is considerably longer 
and -more involved than later methods and need not be con- 
sidered here. Wilsing's method, as originally published, was 
suitable only for orbits of small eccentricity, but Dr. Russell 
has extended it so as to make it applicable to larger 
eccentricities as well. The method of Wilsing and Russell 
is purely analytical, and in time consumed is considerably 
longer than the geometrical methods of Lehmann-Filhes, 
Schwaezschild, and Zurhellen. It consists in finding a 
Fourier's series for the velocity in terms of the elements. It 
should be very useful in special cases, particularly when the 
period is so nearly a year that one part of the velocity-curve 
is not represented by any observations. Another special appli- 
cation suggested by Russell is its use as a method of trig- 
onometrical analysis in disentangling the superposed elements 
of a triple or other multiple system. These longer methods 
will not be further treated here. 

The spectroscopic plates give us a series of irregularly dis- 
tributed radial velocities of the star, and the first problem 
is the determination of the approximate period, which rarely 
offers any difficulty. This should be done as early as possible 
in order to distribute the plates evenly over the revolution 
period and to fill up any gaps which may have been left by 
earlier plates. Different series of plates will differ so widely 
in their arrangement that the determination of each particular 
period is a separate problem. The observation times are most 
conveniently expressed in Julian days and decimals of a day. 
By choosing an epoch, preferably a maximum or minimum 
near the beginning of the series, a probable approximate value 
of the period is easily secured by the use of co-ordinate paper. 
When a number of maxima and minima are well defined, wide 
divergencies in other regions of the curve may generally be 
removed by trying a new period, some sub-multiple of the first 
period assumed. Then by subtracting multiples of the adopted 
period from the observation dates all the observations are 
reduced to a curve covering one revolution beginning at the 
assumed epoch, with the abscissae as times from the adopted 
date of maximum or minimum, and the observed velocities as 
the ordinates. After the preliminary curve is drawn, by noting 
the deviations from the curve at points near the mean of the 
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maximum and minimum velocities, where a change in the 
periodic time will have the greatest effect, advisable changes 
in the preliminary assumed period are readily found by divid- 
ing the deviations of such critical observations, expressed in 
time, by the number of revolutions elapsed. A second curve 
may then be drawn, whose periodic time will generally be very 
close to the true one. For all this preliminary work a calcu- 
lating machine is a great time-saver, though a simple table of 
multiples of the period will be almost as rapid. For stars 
whose periods are almost exact multiples of a day, and short, 
it may be necessary to wait some time before observations can 
be secured over all the velocity curve. The same may be true 
for orbits nearly a year in length, unless the analytical method 
of Wilsing and Russell be used in these cases. 

After the period has been determined as carefully as possible 
the "most probable curve" is drawn, by estimation, through 
the observation points as plotted ; and it is well to prolong 
this curve through a revolution and a half, in order to be able 
to use the ingenious methods of superposition devised by 

SCHWARZSCHILD and ZURHELLEN. 

The errors in drawing this curve have considerable effect 
on the accuracy of the determination of the elements. The 
most probable curve which even a practiced hand may draw 
is not likely to be a perfect representation of the elliptic motion 
which caused it, as it is natural to bend the curve slightly in 
or out at various, points to satisfy the more or less exact 
observations. This difficulty is fundamental in all methods of 
treatment, and for it there seems to be no remedy other than 
by testing rapidly the first orbit found by a trial ephemeris 
and making the needful small changes in the elements, to 
which reference will be made later. 1 



1 Though Zurhellen mentions this unavoidable curve-error earlier in his valu- 
able paper, it appears that he has not duly taken it into account, in the comparisons 
which he makes at the close of his paper, of the values of o> and the eccentricity 
secured by his methods and those obtained by other investigators using other 
methods. His values, it is true, are quite close to the values secured, in some 
cases, by a rigorous least-square solution, but his values have evidently been 
derived from the final curves as published. These, of course, are accurately drawn 
from ephemerides and are very approximately true representations of the line of 
sight component of the velocity in an elliptic orbit. Had he applied his methods 
to the original preliminary curves used by the investigators with whose work he 
makes comparisons, he would doubtless have found values as far removed from 
the final ones as are the preliminary values secured by the method of Lehmann- 
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Figure 2 gives a preliminary velocity-curve as actually 
drawn for the determination of the orbit of the spectroscopic 
binary k Velorum. The observations used were as follows : — 



Julian Day, G. M. T. 


Vel. 


Julian Day, G. M. T. 


Vel. 


2416546.739 


+ 68.5 km 


2417686.59I 


+ 33-8 km 


60.703 


+ 12.9 


9I-57 2 


+ 38.2 


97-65I 


+ 65.7 


92-545 


+ 43-2 


6912.60I 


+ 53-3 


96.480 


+ 46.7 


7587.844 


+ 58.6 


77OI.494 


+ 52-7 


88.788 


+ 57-9 


4I.466 


+ 22.1 


9O.829 


+ 58.5 


46.463 


+ 0.3 


9I.824 


+ 64.8 


49470 


- 7-6 


97788 


+ 65.8 


50.479 


— 8.8 


7609.79O 


4- 62.0 


5I-463 


— 13-3 


54-534 


— 21.0 


53457 


— 19.2 


55-556 


— 19.2 


58451 


— 29.0 


58.570 


— 15-2 


59.460 


— 24.6 


59-545 


— 14-5 







The small circles representing the first four observations, 
which are important in determining the period, owing to their 
distance in time from the later ones, are barred in the diagram. 
The period, U, was assumed to be 116.65 days, and the begin- 
ning of the curve is at Julian day 2416476.0. This is not 
exactly at a minimum, as may be seen from the diagram. 

The conditions of the problem are illustrated in Figure 1. 
The XY-plane is taken as the plane tangent to the celestial 
sphere at the center of motion, the Z-axis being the line of 
sight in which the velocities are measured, and perpendicular 
to the XY-plane. The orientation of the X and Y axes in 
space remains unknown. It will be seen at once that the 



Filhes. Taking the final curve of a 2 Geminorum as published, and using a 
planimeter, I find by the method of Lehmann-Filhes the following values: — 

Lehmann-Filhes. Final Elements, 

e = 0.501 e == 0.503 

w = 2650.36 w = 26 5 °.35 

Using the same procedure for the curve of j3 Auriga as published with Zurhellen's 
paper, I derive e — 0.144, in close agreement with the value from his method 
No. 8, (e = 0.142); for both his method No. 8 and that of Lehmann-Filhes are 
equally sensitive to an erroneous placing of the maximum, which obtains in this 
curve. Shifting the maximum back a distance 0.01 U, as suggested in the paper 
quoted, I find « = o.n (Zurhellen, mean of eight methods, e — 0.1061. 
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problem is somewhat different from that in which the orbit 
of a visual double is determined. In that case the data are 
the observed distances and position-angles at certain epochs, 
which, when plotted, give us a projected ellipse, while the 
spectrographic data are simply the times and the correspond- 
ing velocities directly toward or directly away from the 
observer (-|- when the distance from the Sun is increasing). 
The true velocities or the direction of motion in the orbital 
ellipse, and the size and the position in space of the actual 
orbit, must of course always remain unknown, except in the 
comparatively rare cases where visual and spectrographic 
observations can be combined. It is manifest, moreover, that 
if the orbit-plane were coincident with the XY-plane, such a 
binary could not be detected spectroscopically, as there would 
then exist no component of the velocity in the line of sight. 

The following customary notation will be used: — 
U = Period, 
fi = mean daily motion, 

i = inclination of orbit to XY-plane (not deter- 
minable), 1 
Q, = ascending node on XY-plane (riot determinable), 
= semi-major axis of true ellipse " " 

p = semi-parameter of true ellipse " " 

<<) = angular distance of periastron from ascending 

node, 
u = argument of the latitude, u = v + «>, 
r = radius vector, 
t = the time, 

T = the time of periastron passage, 
v, M, E = true, mean, and eccentric anomalies respectively, 
e = eccentricity, = sin <£, where <j> is the eccentric 
angle, 
V = the velocity of the center of mass of the system 

as a whole (Zurhellen's y). 
f = constant of attraction. 2 

1 While a and i are separately indeterminable, it will be seen later that we can 
find the value of a sin i, from which we may assign a minimum for a. 

2 In the great majority of spectroscopic binaries the spectrum of one component 
only is visible. This does not necessarily mean that the companion is dark, but 
merely that it is considerably fainter. In this case the motion of the bright star 
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By reference to Fig. 1 it will be seen that the distance 
from a point S in the orbit to the XY-plane will be 



r sin 1 sin : 




FIGURE I. 

The spectrograph gives us, however, not the distances from 
the XY-plane, but the velocities of the star's approach or 
recession from this plane, generally expressed in kilometers 
per second, and evidently equal to dz/dt. 

Hence, differentiating, 

dz ... dr , . . du 

-j- = sin 1 sin u — i— + r sin t cos u —. — 
dt dt dt 

From the known laws of elliptic motion we may substitute 

the values 

dr f 

rf7_=y7' sin(M-w) ' 
fy/P f 



du 



= [1 + tfcos (u — <o) ]. 



dt r yp 

These substitutions give for the fundamental equation 



f 



dt 



VP 



= sin i (cos u -\- e cos «). 



(1) 



is determined with reference to the center of gravity of the system and the 
constant of attraction has the value 

km $ 



where k is the Gaussian constant. If both components are visible in the spec- 
trum and the motion of one mass with reference to the other is determined, 
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Now the curve of the observed velocities will have the gen- 
eral appearance of a distorted sine-curve with a maximum and 
a minimum at the points A and B in Fig. 2, but it is evident 
that the curve as a whole is shifted a distance above or below 
the zero axis by an amount equal to the velocity of approach 
or recession of the binary system as a whole, and before the 
observed values can be made purely periodic the velocity 
of the center of mass of the system, which is a constant quan- 
tity for a single system, must be eliminated from the observa- 
tion values. Methods of finding the position of the line repre- 
senting the velocity of the center of mass of the system will 
be noted later, but for a study of some of the properties of 
the curve we will suppose for the present that this velocity 
of the system as a whole is known, and is represented by the 
dotted line in Fig. 2. This will be called the V-axis. 

It is evident from equation ( 1 ) that for a maximum velocity 
cos u = 1, and for a minimum cos u = — 1 ; i. e., at the point 
A of the curve, u = o ; and at B, u = 180 ; or, A and B are 
the points corresponding to the star's passage through the 
ascending and descending nodes of its orbit, respectively. 
Taking, then, A and B as the magnitudes of the curve ordi- 
nates at the points of maximum and minimum reckoned from 
the V-axis, B being taken as a positive quantity, we have 

/ 
A = ^=sin i (1 -\- ecosm), 

f 
B = — 7= sin j'(i — ^ cos w). 

VP 

From these equations it is clear that the term ~y= sin i is the 

half amplitude of the curve in Fig. 2. This will be denoted 
by K, — i. e. 

f „:„v_^_A + B 



v/> 


2 


- -_= sin i e cos o> = 

VP 


A — B 
2 


e cos (o = 


A — B 



A + B 
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The expression for the velocity measured from the zero-axis 

A fj 

may then be written, by placing V = V -\ > 

dz 

— = V -\- K cosu. (2) 

Instead of using the ordinates A and B measured from the 
V-axis to find K, we may of course take the corresponding 

values from the zero-axis, replacing ■ by 



2 



22 

where M t and M 2 are the ordinates of maximum and minimum 
respectively, measured from the zero-axis. 

Given the curve and the above fundamental equations, quite 
a variety of methods may be used to find the elements, all 
requiring about the same amount of time, providing that a 
planimeter is available when using the method of Lehmann- 
Filhes. If the computer is not provided with one of these 
convenient instruments, the methods of Schwarzschild and 
Zurhellen will be more expeditious. Accordingly these 
methods will be treated first. 

schwarzschild's method of finding the periastron 

TIME, T. 

Given M 1 and M 2 , the observed velocities of maximum and 
minimum as measured from the zero axis, draw the line whose 

ordinate is — - — 1 2 . In Fig. 2, M x = + 66.4, M 2 = — 26.8, 

therefore (M 1 + M 2 )/2 = + 19.8. This line will be called 
the mean axis, and is represented in Fig. 2 by the dot-and-dash 
line. Mark on this axis the points U/2 and 3U/2; then lay a 
piece of semi-transparent paper over the curve and copy on this 
the curve with the mean axis, marking also the points o, U/2, 
U, and 3J7/2. Shifting this copy bodily along the mean axis 
for a distance U/2, rotate the copy 180 about the mean axis 
— i. e. turn the copy face downward on the original curve, 
keeping the mean axis in coincidence, and bringing the point 
o or U of the copy over the point U/2 of the curve. The 
curves will then cut one another in general in four points, of 
which two will be the points of periastron and apastron. The 
proper two points will be found without difficulty, for peri- and 
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apastron must be separated in time by one half a revolution, 
and must, moreover, lie on different branches of the velocity- 
curve. To determine which is periastron we have the criteria: 
(1) at periastron the velocity-curve is steeper with regard to 
the axis than at apastron; (2) the curve is for a shorter time 
on that side of the mean axis on which the point of periastron 
lies. By following this procedure we find in Fig. 2 that 
periastron lies at the point -f- 98.8 days from the epoch, giving 
for the dates of periastron J. D. 2416458.15, 6574.80, etc. 

This ingenious method is exceedingly good, though the 
writer prefers the method of Lehmann-Filhes when the 
eccentricity is small. In these cases <*> and T are quite inde- 
terminate, and small errors in drawing the velocity-curve will 
be very troublesome. 

METHOD OF DETERMINING THE ANGLE OF PERIASTRON, to. 

It may easily be verified that V = V + K e cos o> represents 
the position of the mean axis, and accordingly that the ordinate 
of any curve point, referred to the mean axis, may be secured 
from 

1 = 4f — V' = K cos u. 
at 

Now at periastron v = o and hence u = o> ; therefore, if we 
call the ordinate of periastron from the mean axis £ p , we have 

COS to =— • 

A. 

Or, better, we may use the ordinates of both peri- and 
apastron, 



COS *) = 



2K 



Both of these methods are at their best when *> is near 90 . 
From Fig. 2 we take £ p = — 3.2 ; K = 46.6; hence <o = 93°.cj. 
Or, f p = — 3.2 ; U = + 3.9 ; £ p — £ a = — 7.1 ; and o> = 94°4. 
An entirely analogous equation may be used in case the posi- 
tion of the V-axis is known. This is generally most accu- 
rately determined by a planimeter and the method of 
Lehmann-Filhes, to be given later, but may also be found 
from the mean axis by the requirement that the V-axis must 



144 Publications of the 

lie at a distance equal to — ■ K e cos <o above the mean axis. 
Calling the ordinates measured from the V-axis, z, we may 
find 



COS co : 



2K 

Another method of finding o>, due in its simplified form to 
Zurhellen, is best when <o is small. It consists in finding 
the points on the velocity-curve where the true anomaly is 
equal to ± 90°. From the fundamental equation, 



dz A 
ti — 

and 



& = —77 V = K (sin <o + e? cos <o) , for v = — 90 



f 2 = —~ V = K ( — sin <o -\- e cos w) , for v = + 9°° > 

or, t, x = + K sin <o, 

t, 2 = — K sin u>. 

Moreover, at the points v = ± 90 , 

E 1 = — E t , 
M 1 = — M 2 , 
t 1 — T = —(t 2 — T); 

hence these points lie symmetrically in <!; to the mean axis, and 
in time to the periastron. They may accordingly be deter- 
mined by rotating the curve copy 180 about the point of 
intersection of the ordinate of periastron with the mean axis, 
when the two intersections of the copy with the original 
curve will give the points sought. The point v = + 270 may 
be found in like manner. As a check on the location of these 
points draw the lines connecting the points v= -\- 270 and 
v = — 90°, also that from v = — 90 to v = + 90 . These 
lines should cut the mean axis at its intersections with the 
ordinates of apastron and periastron, respectively. Calling the 
ordinates measured from the mean axis ^ for x the point 
v = — 90 , and f, for the point v = + 9°°> then 

sin to = £l ^ 2 , or tan <o = £l ^ 2 
2K C P — U 

This method is at its least advantage in the example given in 
Fig. 2, for here <o is close to 90 . At v ==' — 90 in the figure 
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we find its ordinate from the mean axis to be -f- 46.1, while 
for v = -+- 90 we find — 45.4. 

Hence sin m = 91.5/93.2, and to = ioi°.o. 

Other interesting properties of the velocity-curve and meth- 
ods for finding the periastron-angle and eccentricity are de- 
rived by Zurhellen. For instance, by rotating the curve copy 
about the intersection of the periastron ordinate with the 
V-axis instead of with the mean axis, we may find the curve 
points where the eccentric anomalies are equal to ± 90°. But 
these additional methods will not in general offer any 
advantages over those already given. 

METHODS OF DETERMINING THE ECCENTRICITY. 

I. Perhaps the best is Zurhellen's method of drawing the 
tangents to the curve at the points of periastron and apastron. 
These can be drawn with considerable accuracy, except when 
the periastron falls near a maximum or minimum of the curve, 
for in this case slight changes in the position of periastron 
will cause considerable changes in the inclinations of the 
tangent-lines. The expression for the tangent will be 

dt, 2n dt, 2tt 1: dt, 

dt = U' dM = U 1 — e cos £ ' ~Je' 

From u = v -\- <o, and by introducing the known values 

cos E — sin rf> . cos <j> sin E 

cos v = ~—> sin v — 



1 — e cos E 1 — e cos E 

the original velocity equation will take the form 

dz ,. cos <j> cos <o cos E — sin <a sin E . 

£ = -j- = K cos 4>. = , 

dt 1 — e cos E 

whence 
dt iir ., — cos <4 cos ft> sin E — sin a> cos E -4- e sin 1 

dt U (1 — <?cos£) 3 

At periastron E = o, and at apastron E = ■*■ ; whence 

rf£ — 2ir K cos <j> sin w 



dt v ~ U(i — e) 



•! 



■di, + 2ir K cos <f> sin . 

dt & ~ Uji +e) 2 
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From these values 

tan 

tan 



or ,^-l-I: 



g + l 

In Fig. 2 the value of the tangent at periastron is I04°.8, 
while at apastron it is 6o°.9; whence q = 1.45 and e = 0.18. 

2. Draw the V-axis and determine the ordinates from this 
axis of the maximum and the minimum, and also of peri- and 
apastron. Using z for these ordinates, 

~max ~\~ ^min 

jy v 

~p ^a 

3. The following method is that of Schwarzschild as sim- 
plified by Zurhellen. It is best when <o is small, and to be 
applicable for eccentricities over 0.2 implies the use of 
Schwarzschild's table, of which the necessary portion is 
given below. It depends upon the difference in time of the 
ordinates of the points v = qp 90 , before determined. 

Since tan y 2 E = tan y 2 v tan (45 — y 2 <£), for v = zp 90 , 
we have tan ^= + 1, and therefore E 1 = — (90 — <f>) and 
E 2 — -\- (90 — <f>) . Similarly, 

sin </> sin (90 — <t>) 



.¥, = -(90°-^.) + 
M 2 = +(9o° — <*>) 



sin 1" 

sin <j> sin (90 — <j>) 
sin 1" 



m \/r 3°° f t t \ 00 j sin(i8o° — 2$) 

M 2 — M 1 = ^jj- (t 2 — t t ) = 180 — 2 4> ^ — - — -^> 

U sin 1" 

t 2 — t x may be read off directly from the diagram, and 90 — <f> 
can then be taken directly from Schwarzschild's table. For 
small values of e (<£ < 12 °, e < 0.2) the table may be replaced 
by the shorter formula, 

From Fig. 2 we have t, = 80.0 days, £ 2 =117.8 days, hence 

2 ^ = 0.324; whence from the table, 90 — <f> = 73°.7 and 

e = 0.28. Using the shorter formula, e = 0.27. The di- 
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vergence from the other values which have been secured is 
because <a is so near to 90 , the most unfavorable case for 
this method. 





SCHWARZSCHILD S TABLE FOR 


THE EQUATION 








t. 


! — 'i 






2 V - 


— SIN 2 


ri "^- 2 7T ~ 












~u 






', — '1 


V 


t. — t, 


V 


*•> — ', 


V 


U 


U 


u 


o° 


O.OOOO 


30° 


O.O29O 


6o° 


0.1956 


I 


O.OOOO 


31 


O.O318 


61 


0.2040 


2 


0.0000 


32 


O.O348 


62 


0.2125 


3 


O.OOOO 


33 


O.O380 


63 


0.2213 


4 


O.OOO I 


34 


O.O414 


64 


0.2303 


5 


O.OOOI 


35 


O.O45O 


65 


0.2393 


6 


0.0002 


36 


O.O488 


66 


0.2485 


7 


0.0004 


37 


O.O527 


67 


0.2578 


8 


0.0006 


38 


O.0568 


68 


0.2673 


9 


0.0008 


39 


O.061 1 


69 


0.2769 


10 


O.OOII 


40 


O.0656 


70 


0.2867 


11 


0.0015 


4i 


O.O703 


7i 


0.2966 


12 


0.0020 


42 


O.O751 


72 


0.3065 


13 


0.0025 


43 


O.0802 


73 


0.3166 


14 


0.0031 


44 


O.0855 


74 


0.3268 


15 


0.0038 


45 


O.O9IO 


75 


o.337i 


16 


0.0046 


46 


O.O967 


76 


Q-3475 


17 


0.0055 


47 


O.IO25 


77 


0.3581 


18 


0.0065 


48 


O.I085 


78 


0.3687 


19 


0.0077 


49 


0.1 147 


79 


o.3793 


20 


0.0089 


50 


O.I2I2 


80 


0.3900 


21 


0.0103 


5i 


O.I278 


81 


0.4008 


22 


0.0117 


52 


O.I346 


82 


o.4"7 


23 


0.0133 


53 


O.I416 


83 


0.4226 


24 


0.015 1 


54 


O.I488 


84 


0-4335 


25 


0.0170 


55 


O.I561 


85 


0.4446 


26 


0.0191 


56 


O.1636 


86 


0-4557 


27 


0.0213 


57 


O.I7I3 


87 


0.4667 


28 


0.0237 


58 


O.I792 


88 


0.4778 


29 


0.0262 


59 


O.1873 


89 


0.4889 


30 


0.0290 


60 


O.I956 


90 


0.5000 
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METHOD OF LEHMANN-FILHES. 

It has not seemed advisable to segregate the separate steps 
in this method for inclusion with the other methods given 
before. The method of Lehmann-Filhes will be found very 
generally useful, and little, if any, longer than other methods, 
providing a planimeter is used. 

Given the observations and the velocity-curve, the first step 
is the fixing of the position of the velocity of the center of 
gravity of the system, — the so-called V-axis. This is found 
by the requirement that the integral of dz/dt, — that is, the area 
of the velocity-curve, — must be equal for the portions of the 
curve above the V-axis and below it. By far the easiest 
method of performing this integration is to use a planimeter. 
A line is first drawn by estimation as nearly in the correct 
position to fulfil the requirement of equal areas as can be 
judged by the eye. Placing the planimeter pointer at the 
intersection of this line and the curve, trace out the curve 
lying on one side of the line in, say, a clockwise direction, and, 
on reaching the starting point, traverse the portion of the 
curve on the other side of the line, in the counter-clockwise 
direction. The readings of the planimeter will then give 
directly the difference in the two areas, and the position of 
the V-axis may be shifted slightly to correct the difference. 
A second double circuit by the planimeter will generally be all 
that is necessary to secure an accurate value of the position 
of the V-axis, not depending simply upon the magnitude of 
maximum and minimum as in the previous methods, but sup- 
ported by the entire curve. The increased accuracy is of little 
importance, however, as it will be found best to correct what- 
ever value of V may be adopted by the residuals from the final 
ephemeris. If a planimeter is not available, doubtless the best 
method to equalize the areas above and below the curve will 
be to draw the curve on cross-section paper and make a count 
of the small squares, as suggested by Rambaut. The position 
may also be found by an approximate mechanical integration, 
as advised by Lehmann-Filhes, but in such a case the shorter 
methods explained above for finding the elements would 
generally be given the preference. 

Given the position of the V-axis, the ordinates should next 
be drawn at the points of maximum and minimum velocity 
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(Aa, Bb in Fig. 2). It is at this point that this method is 
weakest and most subject to error, for slight errors in the 
fixing of the positions of these ordinates may easily arise. It 
is best to check the positions of these ordinates by the require- 
ment that the curve portion AaC must be equal in area to CbB, 
while DaA = BbD. 

Considering V as eliminated from the velocities, at the curve 
points C and D, we have dz/dt = o ; hence, from the 
fundamental velocity equation, 

A — B 



cos u = — e cos a> = ■ 



A + B 



The argument of the latitude, u u corresponding to the 
point C, lies on the way from the ascending to the descending 
node, hence sine u ± is positive. Then, 



sin u t ~ 


2\/AB 

" A + B' 






A — B 






A + B 




ile for w 2 at the point 


D 




sin w 2 = — 


2VAB 

A + B 


— — sin u u 


COS M, = — 


A — B 


— COS u, . 



A + B 



Denoting the area of AaC by Z 1; and bBD by Z 2 (it 
will be seen that these areas represent the distances of the 
star from the XY-plane when at the points of its orbit cor- 
responding to the arguments of the latitude u x and w 2 ) and 
the radii vectors by r x and r 2 , we have 

Z x = r-t sin i sin w 1( 

Z 2 = r 2 sin i sin w 2 = — r 2 sin i sin u u 
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Then, since 



r, = 



1 + e cos (u t — w) 
P 



2 1 + e cos (m x + a>) 

1 + e cos u x cos o> — t? sin w x sin «> _ Zj 
1 + e cos w x cos co -|- e sin w t sin &> Z 2 

and, since cos u x = — ^ cos to, 

sin m x — t? sin <o Z, 

sin u 1 -\- e sin <o Z 2 

Whence, 

Z. + Z 

e sin to = -=2 =i- sin m, , 

z, 2 /-j 

which may be written in the form 

2VAB Z 2 + Z, 



e sin to = 



A + B Z 2 — Z x 



As it was found earlier from a consideration of the values 
A and B that 

A — B 



e COS ft) : 



A + B 



these two equations give all that is needed for a determination 
of e and «>. The values of A and B are taken from the curve, 
and the areas Z t and Z 2 are quickly integrated from the curve 
portions AaC and bBD by means of a planimeter (the latter 
area is considered negative in sign and the unit of area used 
is entirely immaterial as only the ratios of the areas are 
required). 

At the time of passing through periastron u = <a, whence 
from the fundamental equation 

= K (1 -f- e) cos to. 



dt p 



Knowing K, e, and m, determine this value, and the abscissa 
of the point on the curve which has the above value as ordinate 
will give the value of T when properly combined with the 
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epoch used for the beginning of the curve. As at the points 
A, B, and A 3 , w equals o, ir, and 2w, respectively, and as o> is 
known, there will be no ambiguity as to the position of u = <o. 
We may instead determine the eccentric anomaly at the 
point C where u x is known, and then find the time from peri- 
astron, t x , corresponding to this value of n x by the well-known 
formulae 

tan £, / t p , M i — w 

L = a/— - tan — , 

2 ^ i+e 2 



2tt 



T = t x 



E x — e sin £, 



The equation for the mean daily motion may be transformed 

as follows : — 

_ 2tt _ / A j- B V7~ I 

U a I 2 sim of 

from which the following equation for the product of the 

separately indeterminable a and sin i may be derived. The 

factor 86,400 must be introduced because the unit of time for 

A and B is the second, while for /x it is the day. 

A _i_ B 

a sin i = 86,400 — — — V 1 — <? = [4.13833] K U Vi — e 2 . 

2/i 

In the use of the method of Lehmann-Filhes slight inaccu- 
racies in drawing the preliminary curve when the eccentricity 
is small are less troublesome than in the methods of super- 
position given before. For eccentricities under 0.10 all meth- 
ods are apt to give too large values of the eccentricity, and it 
is well in such cases to reduce the first value found considerably 
and test a small value first. 

EXAMPLE OF METHOD OF LEHMANN-FILHES (FROM FIG. 2). 

A = 46.3 K = — = 46.6 

2 

6 = 46.9 Z x = + 0.168 
A + B = 93.2 Z 2 == — 0.259 
A — B = — 0.6 Z 2 + Z x = — 0.091 

Z 2 — Z, = — 0.427 
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log VAB 1.6684 
colog(A-j-B) 8.0306 
log (A — B) 9778211 



log Z 2 + Z 1 8.959011 
colog Z 2 — Zj 0.369611 
0.0000 



sin w. 



log 2 O.3OIO 


e sina> 


9.3286 


sin u x 0.0000 


ecosw 


7.88o8n 


cos u x 7.8088 


tan <o 


i.4478n 


Ml 8 9 °.6 3 


«) 


92°.04 




\oge 


9.3289 




e 


0.21 


log (1 + e) 0.0828 


const, log 


4-I383 


cos<o 8.55 1 911 


K 


1.6684 


log if 1.6684 


U 


2.0669 


0.303 m 


Vi—e 2 


9.9902 


ordinate p = — 2.o km 


a sin i 


7.8638 


t p = 98.4 days 


. =73, 


ooo,ooo km 



= — cos u x ) 



T = J. D. 2416457.75 

(The longer method gives 6458.0.) 

IMPROVED ELEMENTS. 

Whatever the method used to determine the elements, it will 
always be necessary to test the elements as derived and make 
small changes in them before they can be regarded as a satis- 
factory representation of the observations. To do this a test 
ephemeris is rapidly computed for ten or twelve points of the 
curve. V being, as before, the velocity of the system as a 

A B 

whole, and placing V = V -\- K e cos m = V -\ , the 

2 

necessary formulae for the ephemeris are 

dz 



-^=V> + Kcosu, 



tan 



(u- 



o=Vj 



+ e 



tan 



E, 



M = E — e sin E, 

M = n(t—T). 

Collecting the elements as derived by the method of 
Lehmann-Filhes, with the necessary constants, we have: — 
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U — 116.65 days ii° = 3°.o86i 



,V4± 



e 



e — o.2\ log\ — - — =0.0926 

T = 2416458.0 = — 0.3 

K = 46.6 V' = + 20.4 

(1) = Q2°.04 

F = -f 2o.7 km 

Then, at the point where t = 20 days, on Fig. 2, we have 
the following sample from the test ephemeris: — 

t = 2416496.0 
t— 7=38.0 y 2 (« — o.) =79°.23 

M =148°. 13 M 0) =158.46 

£ = 153.50 (From Astrand's Hiilfstafeln) W = 250.50 

tan y 2 E = 0.6281 nat cos u = — 0.334 

tan y 2 (u — «)) = 0.7207 F' + /C cos u = cte/* = + 4.8 km 

In this way a test ephemeris was rapidly computed for the 
points o, 10, 20 days, etc., of Fig. 2. The resulting curve, to 
save confusion, is not drawn in, but its position is represented 
by the heavy black dots. A study of the discrepancies shows 
that e is a little too large, the branch on the apastron side 
being somewhat too sloping, while the other slope is a little 
too steep. When e is changed we shall also need to make 
small changes in T and o>. Making one or two changes in 
this way, and testing at critical points of the curve, the fol- 
lowing elements were finally decided upon as best satisfying 
the observations : — 

V = + 2i.9 km 

U = 116.65 days 

e = 0.19 

K = 46.5 
u) = 96°.23 
T = J. D. 2416459.0 
a sin i = 73,ooo,ooo km 

The correction to the value of V was found last of all from 
the residuals of the final ephemeris by the simple formula 
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, where n is the number of observations and v the resid- 

n 

ual, o — c. The residuals from the final ephemeris and the 
final curve may be found in Lick Observatory Bulletin, No. 
122, 1907. 

Combining the best results from the methods of Schwarz- 
schild and Zurhellen, we would have the following 
elements : — 

7 = 2416458.15 
<o = 94.4 
e = 0.18 
K = 46.6 
V' = + 19.8 



] °syl~zv=°- 



0790 

A test curve from these elements was also computed, and its 
position is represented by the small crosses at the points 5, 15, 
25, etc., days. These elements are seen also to need small 
corrections, but on the whole are rather nearer the observa- 
tion curve than are the values secured by the method of 
Lehmann-Filhes, as the position of periastron is about the 
most favorable for the determination of the eccentricity by 
Zurhellen's methods of tangents. 

With some experience it is possible to secure elements which 
would be little, if any, bettered by a least-square solution, 
unless a large number of well-determined velocities are avail- 
able, extending over a considerable interval of time. For 
completeness, however, the least-squares formulae are given 
here as derived by Lehmann-Filhes: — 

dz 
8 — = (cos u -\- e cos <o) 8A' 

. „ r sin u sin v , . . , . 

-\- K [cosw — (2 + ^cosw)] tie 



K (sin u -\- e sino>) 8 



O) 



— if sin u (1 + e cos^) 2 . (t — T) 



(1— «■)* 



+ sinw (1 + <?cos v) 2 ST. 

(i—e 2 y 
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In case both spectra are visible and the motion of one star 
is determined with reference to the other / = k yjm -\- m', 
and an equation for the relation of the masses is found from 
the relation 



giving 



k y ' m + m' 

/*= 5 : 

a * 



a 2 (a sin i) 1 

m -\- m = 



k 2 sin 3 i 



where A, B, and a sin i are to be expressed in terms of the 
mean distance of the earth from the Sun. If the velocity of 
the center of mass of such a system has also been determined, 
then from this the relation of both masses is known, for each 
mass may then be represented separately in the form 

constant 



sin° 1 



Further data as to methods of treating the comparatively 
rare cases where both spectra are measurable are given in 
Zurhellen's last paper. 

In the still rarer cases where both visual and spectrographic 
observations of a system are at hand, valuable results as to the 
parallax and relative masses of such systems may be secured. 
Reference may be made here to Professor Hussey's deter- 
mination of the parallax of 8 Equulei {Lick Observatory Bul- 
letin, No. 32), Professor Wright's value of the parallax of' 
aCentauri (Lick Observatory Bulletin, No. 60), and to Dr. 
Campbell's discussion of the variable radial velocity of Sirius 
(Lick Observatory Bulletin, No. 70). 

The D. O. Mills Expedition, 

Santiago, Chile, October, 1907. 



